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Abstract 

The analytic proof of mode stability of the Kerr black hole was provided by 
Whiting. In his proof, the construction of a conserved quantity for unstable mode 
was crucial. We extend the method of the analysis for the Kerr-de Sitter geometry. 
The perturbation equations of massless fields in the Kerr-de Sitter geometry can 
be transformed into Heun's equations which have four regular singularities. In this 
paper we investigate differential and integral transformations of solutions of the 
equations. Using those we construct a conserved quantity for unstable modes in the 
Kerr-de Sitter geometry and discuss its property. 
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1 Introduction 



One of the most non-trivial aspects of perturbation equations for the Kerr-de Sitter 
geometry is the separability of the radial and angular parts. Carter first showed that the 
perturbation equation for a scalar field is separable in the Kerr-Newman-de Sitter geom- 
etry [|J. This observation was extended for spin 1/2, electromagnetic fields, gravitational 
perturbations and gravitinos for the Kerr geometries and even for the Kerr-de Sitter class 
of geometries. These perturbation equations are called Teukolsky equations ||. Except 
for electromagnetic and gravitational perturbations, the separability persist even for the 
Kerr-Newman-de Sitter solutions. 

An important application of the separability is the proof of the stability of the black 
hole. The proof of mode stability of the Kerr black hole was provided by Whiting P| in 
1989. The proof is more complicated than one for the Schwarzschild black hole, because 
in the Kerr geometry there is no Killing vector which is timelike everywhere in the region 
exterior of the outer horizon. In his proof, he skillfully used differential and integral 
transformations of solutions of perturbation equations, and obtained a conserved quantity 
which is well-defined for unstable modes which are purely incoming on the outer horizon 
and purely outgoing at the infinity, and are characterized by positive imaginary part 
of their frequencies. Then he showed that the positivity of the quantity bounds the 
magnitudes of the time derivative of perturbations. 

On the other hand, it was shown that the perturbation equations in the Kerr geometry 
are obtained from those in the Kerr-de Sitter geometry in the confluent limit. In other 
words, a irregular singularity at infinity of the equation in the Kerr case is separated into 
two regular singularities of the equation in the Kerr-de Sitter geometry by cosmological 
constant. In a series of paper j|, ||, |J, Suzuki, Takasugi and the author have constructed 
analytic solutions of the perturbation equations of massless fields in the Kerr-de Sitter 
geometries. We found transformations such that both the angular and the radial equations 
are reduced to Heun's equation ||. The solution of Heun's equation (Heun's function) 
is expressed in the form of a series of hypergeometric functions, and its coefficients are 
determined by three-term recurrence relations 0. The solution of the radial equation 
which is valid in the entire physical region is obtained by matching two solutions which 
have different convergence regions in the region where both solutions are convergent. We 
examined properties of the solution in detail and, in particular, analytically showed [|| 
that our solution satisfies the Teukolsky-Starobinsky identities ||. There are similarities 
between the procedures for solving the perturbation equations in the Kerr and Kerr-de 
Sitter geometries. Thus differential and integral transformations of Heun's function may 
be useful for studying mode stability of the Kerr-de Sitter geometry. 

In this paper, we investigate differential and integral transformations of solutions of 
massless perturbation equations in the Kerr-de Sitter geometry. These transformations 
map a solution of a Heun's equation to a solution of another Heun's equation which has dif- 
ferent parameters from original ones. Differential transformations include the Teukolsky- 
Starobinsky identities as special cases, and also other transformations for angular func- 
tions. We will consider two specific integral transformations because we do not know 
systematic way to study integral transformations of Heun's function. Then we will apply 
the integral transformation to a solution of the radial equation. It will turn out to be 
possible to make the transformations of radial functions only for unstable modes which 
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are purely incoming on the outer horizon and purely outgoing on the de Sitter horizon, 
and have positive imaginary part of their frequencies. A conserved energy integral will be 
constructed from the transforms of the angular and radial functions for unstable modes, 
and we discuss properties of the quantity. In the Kerr limit, these transformations and 
the conserved quantity are coincide with those given in Ref. by Whiting in order to 
prove mode stability of the Kerr black hole. 



2 The Teukolsky equations for the Kerr-de Sitter 
geometry 

We consider perturbation equations for massless fields in the Kerr-de Sitter geome- 
tries. In the Boyer-Lindquist coordinates, the Kerr-de Sitter metric has the form, 



(l + a)V 

where a = Aa 2 /3, p 2 = pp*, p = r + ia cos 9 and 



dr 2 d9 2 \ A 9 sin 2 # 

a7 + a^J ~ (i+cov 1 

(dt- a sin 2 9dtp) 2 , (2.1) 



ds 2 = -p 2 f ^ + ^ i - ,r " i adt - (r 2 + * 2 )m 2 

A 



A r = (r 2 + a 2 ) ( 1 - -^r 2 ) - 2Mr 



a 



a 2 
ct 

= — ( r _ r+ \( r _ r _\l r — r ' V r — r_) , 
a 2 

A e = l + acos 2 9. (2.2) 

Here A is the cosmo logical constant, M is the mass of the black hole, and aM is its 
angular momentum. 

We assume that the coordinate dependences of the perturbation of the field have form 
$ s = e -H w *- m w R s (r)S s (9). Then angular Teukolsky equation for massless field with spin 
weight s is given by 



dx dx a 

1 + a 

- 2s(am — (1 + a)^)x — - — 

a 



+ 



2s(am - (1 + a)£)x - - + ^ £ 2 + 2m(l + a)^ + s 2 



1 + ax 2 

(1 + a) 2 m 2 (1 + a)(s 2 + 2smx) 

(1 + ax 2 )(l — x 2 ) 1 — x 2 



>S a (x) = 0, (2.3) 



where x = cos 9 and £ = auj. The separation constant \ s is an even function of s as shown 
in Ref. [[§] . This equation has five regular singularities at ±1, ±^7= and 00. We define the 
variable z by 



2 x-^r 



(2.4) 



2 



Then the singularities are transformed to z = 0, l,z s ,Zoo, and oo where z s - 

and = ~ '^2^? • The singularity at z = z^ which corresponds to x = oo can be 
factored out by the transformation 



S s (z) = z c \z - \f\z - z s f*(z - z^fsiz), 



(2.5) 



where C\ = 5 x (m-s)/2, C 2 = 5 2 (m + s)/2 and C 3 = 5 3 | (yf £ - v 7 "" 1 _ * s )> (^1^2, <$3 
±1). Now fs(z) satisfies the equation 



dz 2 



+ 



2Ci + l | 2C 2 + 1 | 2C 3 + 1 



z - 1 



z - z s 



d 

— + 



P+P-Z + u 



dz z(z — l)(z — z s ) 



fs(z) = 0, (2.6) 



where 



p± =c 1 + c 2 + c 3 ±c; + i 

—i ' 



u - 



4^ 



A s - 2i^a~ + 2(1 + a)(m + s)£ - (1 + iy/a) 2 (2CiC 2 + d + C 2 ) 

2 



- 4i>/a(2C r 1 C 3 + Ci + C 3 ) - ^- [4a + (1 + i^) 2 

s 2 1 
+ — (1 — iy/a) 2 + 2imsy/a(l + «v^) r- 



(2.7) 
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Equation (2.6) is called the Heun's equation which has four regular singularities. 
Next we consider the radial Teukolsky equation which is given by 



d 



A- S —Ai 



s+l 



d 1 

+ 



r dr r dr A, 



1 + a) 2 K 2 - i a (l + a)K 



dA r 
dr 



2a 



+ 4zs(l + a)ur - —(s + l)(2s + l)r 2 + s(l - a) - X s }R s (r) = 0, (2.9) 



with K = u(r 2 + a?) — am. This equation has five regular singularities at r±,r'±, and 
00 which are assigned such that r± — > M ± -\/M 2 — a 2 — Q 2 = r± and — > ±a/ \/a in 
the limit a — > 0(A — > 0). And the coefficients of the equation are complex for spin fields 
(s ^ 0). We assume that the cosmological constant is sufficiently small so that all r±,r± 
are real. By using the new variable 



(2.10) 



( r + 


-r'\ 


/ r 


— r_ \ 




— r_ / 


V r 


-rl / 



Eq. (|2.9|) becomes an equation which has regular singularities at 0, 1, z r , z^ and 00, 



( r + 


-r'\ 


(r' + 


— r_\ 


V r + 


— r^J 




- r'_) 



r _ r 



r 1 — r_ 



To proceed further, we define the parameters 



A± = -{-s±(2 ai + s)}. (^=1,2,3,4) 



(2.11) 



Here a's are purely imaginary numbers defined by 

a 2 (l + a) (uo(r\ + a 2 ) — arnj 



a i 



«2 



a 3 



a i 



a(r' + — r + )(r'_ — r + )(r_ — r + ) ' 

a 2 (l + a) (cu(r 2 + a 2 ) — am) 
i — 

a(r+ — r_)(r / _ — r_)(r + — r_) ' 

a 2 (l + a) (co(r' + 2 + a 2 ) — am) 
I — 

a(r_ — r' + ){r'_ — r' + ){r + — r' + ) ' 

a 2 (l + a) (u{r' 2 + a 2 ) - am) 
i — 

a(r_ — r'_)(r' + — r / _)(r + — ri.) ' 



(2.12) 



and the relation ai + a 2 + a 3 + a 4 = is satisfied. Again we can factor out the singularity 
at z = by the transformation as 

R s (z) = z D \z - \) D \z - z r ) D *(z - Zoo ) 2s+1 f R (z), (2.13) 

where Di (i = 1,2,3) is either D i+ or D^. Then, fii(z) satisfies the Heun's equation as 



dz 2 



+ 



7 5 e 
+ + 



d 



o + o_z + V 



dz z(z — l)(z — z r ) 



f R (z) = 0, 



^z z — 1 z — z r 
where 

7 = 2D 2 + s + l, 6 = 2D 1 + s + l, e = 2D 3 + s + l, 

a± = D 1 + D 2 + D 3 + D^ + 2s + l. 

The parameters 7, 5, e and a± satisfy the following relation, 

7 + 5 — 1 = a + + cr_ — e, 



(2.14) 



(2.15) 



(2.16) 



which is required for Eq. (|2.14j ) to be a Heun's equation. The remaining parameter v is 
given by 

2a 4 (l + a) 2 (r+ - r'f(r + - r'_) 2 (r_ - r'_)(r' - r'_) 

v = — — 

a 2 T> r + — r_ 

x I — cj 2 ri(r + r_ — 2r + r' + + r_r^_) + 2auj(auj — m)r-(r+r' + — r 2 ) 
—a 2 (auj — m) 2 (2r_ — r + — r' + ) 



2isa 2 (l + a) 



uj{r_r'_ + a 2 ) — am 



(r + — r_)(r+ — r / _)(r_ — r'_) 



2r 



/ 2 



+(s + l)(2s + l) 

(r + — r_)(r+ — r_) 

-2D 2 {z r Dx + L> 3 ) - (s + 1) [(1 + 2 r )L> 2 + Z r £>l + £3] 



Here T> is the discriminant of A r = 0, 



V 



(r + — r__) 2 (r + — r' + ) 2 {r + — r'_) 2 (r_ — r' + ) 2 {r_ — r'_) 2 {r' + — r'_) 2 
16a 10 ' 



'I -a) 



a- 



M 



1 — ftW 



+- 



ft 



-27M 4 + 36(1 - ft)MV - 8(1 - ft) V 



16ft 2 



It should be noted that the equation fl2.9p for s = is the perturbation equation for 
a conformal scalar field which satisfies 0<f> = \R(j). In the case of an ordinary scalar field 
which satisfies □</> = 0, the term — ^fr 2 in Eq.( [2.9| ) is absent. 



3 Differential transformations 

We consider a differential transformation of a Heun's function f(z) which satisfies 

r 9 2 



M 2 (7,5,e;ft,/3;g) /(z) = < z(z - l)(z - a H ) 



dz 2 







+ [y(z - l)(z - a H ) + 8z(z - a H ) + ez(z - 1)] — + ft^ + q jf(z) 
Differentiating this equation N times, we obtain 

N f 



0. (3.1) 



j^j M z (^,8,e;a,ft;q) f(z) = ^z(z - l)(z - 
+ [( 7 + iV)(z - 1)(2 - a H ) + (8 + N)z{z - a H ) + (e + N)z{z - 1)] 



dz 



+ 



3iV + (2ft + 2ft - 1)N + aft\ z + q - N(N - l)(a H + 1) - N((a H + 1) 7 + a H 8 + e) 
+N(N-l + a)(N-l + ( 3)(4-) 1 )(4z) /(*)• (3-2) 



Therefore 



JV 



.dz) ^ Z 



(3.3) 



with N = 1 — ft or 1—ft formally satisfies another Heun's equation in which the parameters 
(7, 8, ■ ■ •) are replaced by 7 = 7 + iV, 8 = 8 + N,e = e + N, q = q- N(N - l)(a H + 1) - 
N((a,H + 1)7 + Off<5 + e), and ft, ft which are determined by relations a+ft=a+ft+ 3N 
and aft = 3N 2 + (2a + 2ft - 1)N + aft. Here we choose N = 1 - a. Properly N should 
be an positive integer. In both angular and radial cases, we can take parameters so that 
N = |2s|. These transformations are the Teukolsky-Starobinsky relations |§. In Ref.|| 
it was shown that the radial Teukolsky-Starobinsky relations are satisfied by using an 
analytic solution expressed by the form of a series of hypergeometric functions. 

In the angular case, N = \s — m\ and iV = \s + m| are also possible. We now 
explain the case with 8\ = 82 = —83 in the definitions of Cj (i — 1, 2, 3), explicitly. Then 
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N = Si(s — m) and thus we can take 5i(s — m) = \s — m\. The transform of angular 
function fs(z) = fs(z) satisfies the Heun's equation with the following parameters, 



7 = 1, 5 = 25^ + 1 
p+ = 1 + \s — m\, 



Si 



i ~t_ £ + (1 — iy/a)m 



+ 1, 



-<5iz 



1 + a 



a 



(3.4) 



u = u — \s — m\ 



z s (5 1 (m + s) + 1) — <M 



1 + a 



£ — \/«™ — is ] + 1 



We next define S s (z) from /s(z) through thejnverse procedures to those used to derive 
fs(z) from S s (z). By setting 7 = 2C\ + 1,5 = 2C2 + 1 and e = 2C3 + 1, we obtain 
C*i = 0, C*2 = S\s and C3 = — y + (1 — iy/a)m . Using these (7j as the exponents 

at the singularities, S s (z) is given by 

S s (z) = z^(z - lf 2 (z - z s f 3 (z - Zoo )f s (z). (3.5) 

This new angular function S s satisfies a similar equation to the angular Teukolsky equation 
( |2.3|) satisfied by S s ; 



~ x2 )i l + ax 2 )4- + A s - 2ax 2 - (1 + a) 2 ^ — ^ 2 + 2(1 + a) 2 -^ — % 
ax ax 1 + ax z 1 + ax z 



-ail + a)- 5-m 

1 + ax 1 



2 fl + a )^s 2 \sjx) = 0, 



1 — x 



(3.6) 



where we used (|2.4|). We note that the operator acting on S s in the above equation is 
invariant under s — > —s because A, is an even function of s. 



4 Integral transformations 

We construct integral transformations of Heun's function f(z) which satisfies Eq. ( |3.1| ). 
The integral transformation which maps a solution of a Heun's equation to a solution of 
another Heun's equation is given by 



f(z) = [dt f~\t - l) 5 - 1 ^ - a H y~ l K{z,t)f{t). (type A) 



(4.1) 



The function f(z) is a solution of M z (Jy, 8, e; oc, ft, q) f(z) = if the kernel satisfy the 
condition 

M 2 (7, 1 e; a, ft q) - M t ( 7 , 8, e; a, ft q)] K(z, t) = 0, (4.2) 
and the surface term of the integral 

f(t - l)\t - a H y \k(z, t)^-- ) = W(z, t), (4.3) 



dt 



dt 



vanishes at the ends of C. 

In this paper we consider the kernels which depend on z and t only by the form of 
their product zt = ( thus K(z,t) = K((). Then Eq.(4.2) can be rewritten as 

d 2 



t 



C(C-a H )^+[(a + (3 + l)C-a H ^)— + a(3 



d 



d 2 d 
C(C - «^)37? + ((" + £ + 1)C - o-Hl) -T7 + a(3 



d 



;i + a H )(7 - 7) + - 5) + e - ejC^ + g - g WC) = 0. 



In particular, we require the kernel to satisfy the following equations, 

d 2 



(4.4) 



C((-a H )^ + {(a + f3 + l)(-a H ^) — + a(3 



d 

dC 



d 2 d 
C(C -a H )— + ((a + P + 1)C - a H7 ) j7 + «/3 



/f (0 = o, 
#(0 = o, 



(4.5) 



(1 + a H )(7 - 7) + a^(5 - 5) + e - e = 0, 
9 = ?■ 

There are four sets of solutions of these equations. Using invariance of the equations 

under a < — > (3 and a < — > (3, we can take the solutions without loss of generality as 

K(C) = (C-a H y , (4.6) 

5 = 7, 13 = f3, 7 = a, 5 = /3-e + l, e = /3-<J + l, g = g. (4.7) 

Next we define another type of integral transformation, 



/(*) 



a a 



dt t^\t - l) s -\t - a H y- l K{z,t)f{t), (type B) 



(4.8) 



where it should be noted that the integral region depends on z. In order for f(z) to satisfy 
Heun's equation, the following condition must hold instead of the conditions 



W(z,t)\f + z(z-l)(z-a H )l-^- 



^LK(z 1 t)t^ 1 (t-l) 5 - 1 (t-a H y 1 f(t) 



a H dK(z, t) t7 _ 1(t _ 1)5 _ 1{t _ aH)e _ ! 



z 2 dz 

+ [j(z - l)(z - a H ) + 5z(z - a H ) + ez(z - l)p±K(z,t)fi-\t - \Y~\t - a^) e ~V(t) 

z 

= o. 



(4.9) 



We can show by using the kernel in Eq.( [4.6| ) that f(z) is a solution of the Heun's equation 
with the parameters in Eq. (|4.7|) , provided that this condition is satisfied. 
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5 Conserved energy integral for unstable modes 



In this section we are going to construct a conserved quantity for unstable modes 
in the Kerr-de Sitter geometry by using the differential and integral transformations of 
Heun's function given in section 3 and 4. We will assume that the cosmological constant 
is sufficiently small so that r'_ <C r_ < r + <C r' + . 

First, it can be easily checked that the differential transformations of the angular 
function for the Kerr-de Sitter black hole which we provided in section 3 coincide with 
those given by Whiting in the Kerr limit A — > (a — > 0). Indeed since the Teukolsky 
equations for the Kerr black hole have the forms of confluent Heun's equation 0], this 
differential transformation in the Kerr limit becomes one for the confluent Heun's function. 

Next we consider the integral transformations of the radial function fn(z) which satis- 
fies Eq. ([2.14j) . We choose D\ = —ai — s, D 2 = — a 2 — s and D 3 = —a 3 — s as the parameters 



included in the radial equation ( 2.14 ). Then we have 
a + = —2s + 1, cr_ = 2a 4 — s + 1, 

7 = -2a 2 -s + l, <5 = -2ai-s + l, e = -2a 3 - s + 1. (5.1) 

The unstable modes are purely incoming on the outer horizon z = 1 and purely outgoing 
on the de Sitter horizon z = z r , 

RJz) ~ (z-l)— ai , (z~l) 



z 



a:', 



[z ~ Zr) (5.2) 



f R {z) ~ (z-l)°, (z 



z \ 2a 3 +s 

1-- , (z~z r ) (5.3) 



or equivalently 



z r 

where we used the parameters determined above. Furthermore the unstable modes are 
characterized by having positive imaginary part of frequency u. In our discussions below, 
we choose the integral region as C = (1, z r ) in the integral transformation of type A. 
We first examine the Kerr limit. In the Kerr limit, we find 

«4 — > 7r^ ; ( 5 - 4 ) 



2V^K-r°)' 4 2v^ 



where r° = M ± \JM 2 — a 2 . Then both integral regions of the transformations of type A 
and B become (1, oo) and the kernel (|4.6|) becomes of the Laplace type, 



K(C) ~ e 2 ^ (r +" r - K . (5.5) 

Thus both integral transformations coincide with those used in the proof of mode stability 
of the Kerr black hole in this limit. 

We next consider the boundary term W(z,t) . In type A case, W(z,t) behaves as 

W(z,t) ~ (t-l)- 2ai ~ s+ \ (t~l) 

~ (t-z r )- 2a3 - s+1 ^(t-z r ) 2a3+s . (t~Zr) (5.6) 
at 
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The real part of ai is negative if the imaginary part of u is positive. Thus if s < 0, 
W(z,t) vanishes at t = 1 for unstable modes. On the other hand, W(z,t) does not 
vanish at t = z r . Therefore this transformation is not appropriate. Although there is 
also a possibility that our choice of the integral region C is wrong, we do not adopt the 
transformation of type A in our discussions below. In type B case, W(z,t) behaves as 



W(z,t) ~ (t-1) 



z 



-2ai— s+l 
-2a4+s— 1 



(t~l) 



Z r 

t ~ — 

z 



(5.7) 



Since the real part of is negative for unstable modes, it is clear that the condition (|4.9| ) 
holds for unstable modes from the form of the kernel K(z, t) = (zt — z r )~ 2a4+s ~ 1 . Hence it 
is possible to make the integral transformation of type B for radial functions of unstable 
modes. 

The radial function fn{z) given by the integral transformation of type B of fn(z), 



Uz) 



dt t 



-2a,2 — si 



(t-l)- 2ai - s (t-z r ) 



-2a 3 - 



} (zt — Z r ) 



-2a4+s— 1 



/*(*), 



(5i 



satisfies the Heun's equation with parameters ( |4.7|) , 



a + = —2a2 — s + l, (J- — 2a4 — s + l, 

7 = -2s + 1 = 2D 2 + 1, 5 = -2oi - 2a 2 + 1 = 2D X + 1, 

e = -2a 2 - 2a 3 + 1 = 2D 3 + 1, v = v. 

We define a new radial function R s (z) from f R by 

R s (z) = z b \z - \f\z - z r f\z - z^)j R {z). 
This function satisfies the following equation 



(5.9) 



(5.10) 



^-A r ^- - ^ - s 2 F s {r) - m 2 F m {r) + u?F u (r) + m^F mw (r) - A s | R s {r) = 0, (5.11) 



where 
F s (r) 

F m {r) 
FJr) 



~~?( r + ~ r-)(r'_ - r_)- — — + ^(r+ + r_V 
a z r — r_ cr 



4a 4 (l + a) 



r — r_ 



a 2 {r + — r_) 2 (r / f — r_) 2 (r'_ — r_) 2 (r — r+)(r — r+)(r — rL) 

(r-r+)}, 

a 2 (l + a) 2 (r — r_) 



x<{ a(r + - r_)(r + + r_) 2 + (1 - a)a 2 - 2ar 2 



a 3 (r + — r_) 2 (r+ — r_) 2 (r'_ — r_) 2 (r — r + )(r — r+)(r — rL) 
x <{ (r 2 - r 2 ) 2 [a 2 (l + a) - a(r + + r_) 2 

2(1 - a)a 6 + a 4 ((l - 4a + 5a 2 )r 2 - 4a(l - 2a)r+r_ + (1 - 8a + 5a 2 )r 



+2(r H 



2aa 2 (1 - a)r% + 2(1 - a)r\r_ + (3 - 7a)r 2 + r 2 _ - 8ar + ri + (1 - ha)r 



+ 



+cr(r+ + 4r+r_ + 9r\_r 2 _ + 8r^_ri + 5r 2 rl + r_) (r - r + ) 

4a(l - a)a 6 + a 4 ((l - a) 2 r 2 - 2(1 - a) 2 r + r_ + (1 + 6a - 15a 2 )r 2 ) 

-2cra 2 ((l - a)r% - 4(1 - a)r+ri + (1 + 7a)rt) 

+a 2 (r ( \_ + 2r\r_ + 3r 4 r 2 - Ar 3 + r 3 _ - 5r 2 r 4 . - Qr + r 5 _ + r 6 ) (r - r + ) 2 

4a 3 (l + a) 2 (r-r_) 
a 3 (r + — r_) 2 (r' + — r_) 2 (r'_ — r_) 2 (r — r + )(r — r+)(r — r'_) 

x|a(r 2 -r 2 ) 2 -(1 + a)a 2 + a(r+ - r_) 2 

r_) [ - (1 - a)a 4 - 4a 2 (r + + r_) V 
+cra 2 ((l - 3a)(r+ + r_) 2 + 2(1 + a)r*)] (r - r+) 

-2a(a 2 + r 2 ) [a 2 (l - a) - 2ar 2 1 (r - r + ) 2 1, 



(5.12) 



Here we used Eq.( |2.10| ) for rewriting the equation in terms of r. All the coefficients of 
this equation are real for any spin weight s in contrast with the radial Teukolsky equation 
791). We note that RJr) behaves near r ~ r + and r\ as 



(5.13) 
(5.14) 

dt 2 
d 2 



R s (r) ~ (r-r + )- ai ~ a \ (r ~ r+) 

~ (r-r;)~ ai - a4 . (r~r+) 

We construct the function $ s from S s (^) and R s (r) as 

$ s = e -^- m ^^(r)^(0), 

then this function satisfies 



d d 1 d 2 ^ 

dr r dr ~^ sin 9 d9 S ' ( + a os qq 



2„2 



+ 



F TO (r) + a(l + a 
F s (r) + (l + a 



( (1 -cosfl) 2 
1 + a cos 2 

1 + cos 9 
1 + a cos 2 # 



+ 



F u (r) - (1 + afa 
F muj + 2(1 + a) 2 a 



sin 2 fl 



1 + a cos 2 
1 — cos 9 



2a — + cos 
V a" 1 



1 + a cos 2 # 
= 0. 

-s. 



dtdip 
(5.15) 



All the coefficient in this equation are real and invariant under s — 

Finally we obtain a conserved energy integral from this equation in the following form: 



drd9dip sin^j 
+ (1 + a cos 2 



Fu{r) - (1 + 



\2 2 

a) a 



sin 2 9 



1 + a cos 2 9 



<9$ s 


2 r- 


89 


+ 



+s' 



F s (r) + (1 + a) 



F m (r) + a(l + a 

1 + cos 9 
1 — cos 9 



(1 



dt 



COS I 



9$. 



1 + a cos 2 6* 

/ y2 



9$ 



<9r 

2 



dip 



+ 2a I — + cos 2 6> 



(5.16) 
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where r integration is performed over (r + , r' + ). From Eq.( |5.13 ), it can be understood that 
this integration is finite for unstable modes, provided that the cosmological constant is 
sufficiently small. 



6 Summary and discussions 

Solutions of the perturbation equations of massless fields in the Kerr-de Sitter ge- 
ometries can be obtained by using Heun's functions. In this paper, we constructed the 
differential and integral transformations of Heun's function. The differential transforma- 
tions of the angular and radial functions include the Teukolsky-Starobinsky identities as 
special cases. Although we don't know generic way to study integral transformation of 
Heun's function, we provided two types of integral transformations. In the Kerr limit, 
these differential and integral transformations coincide with those considered in the case 
of perturbations for the Kerr geometries f3j. From the equation satisfied by the trans- 
form of the perturbation, we have succeeded in obtaining a conserved energy integral for 
unstable modes. 

In the proof of mode stability of the Kerr black hole, it is crucial that the conserved 
quantity which is obtained by similar procedures to those used in this paper is positive 
definite. From the positivity, it is concluded that the value of the conserved quantity 

bounds the integral of ^f- and thus unstable modes cannot exist. However the con- 
served quantity obtained for the Kerr-de Sitter black hole is not positive definite for the 
sufficiently small ( but non-vanishing ) cosmological constant because the coefficient of 

in the quantity, which vanishes in the Kerr limit A — > and the Schwarzschild-de 
Sitter limit a — > 0, can become negative in the situation considered here. Therefore 
we cannot rule out possibility that there are unstable modes. We also point out that al- 
though is globally null in the metric derived from the equation in the Kerr geometry || 



which can be obtained from Eq. (|5.15| ) in a — ► limit, it does not hold in the metric 



derived from Eq.( |5.15[) with non- vanishing cosmological constant. However we think that 
the procedures performed here are natural extensions of those used in the Kerr case. It 
may be possible to improve the analysis which we provided here. To the end, we think 
that the systematic study of integral transformations of Heun's function will be required. 

The analyses given in our previous papers j|, [5|, || and here are applicable to the 
case of the Kerr-anti-de Sitter geometry similarly. We hope that those analyses may give 
deeper insight for the correspondence between quantum gravity in anti-de Sitter space 
and conformal field theory defined on the boundary jB|. 
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